Thomas-Fermi scaling in the energy spectra of atomic ions 
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The energy spectra of atomic ions are re-examined from the point of view of Thomas-Fermi 
scaling relations. For the first ionization potential, which sets the energy scale for the true discrete 
spectrum, Thomas-Fermi theory predicts the following relation: Eioniz = Z'^ N~'^^^ g{N / Z), where 
Z is the nuclear charge, is the number of electrons, and <; is a function of N/Z. This relation does 
not hold for neutral atoms, but works extremely well in the cationic domain, Z > N . We provide 
an analytic expression for g, with two adjustable parameters, which fits the available experimental 
data for more than 380 ions. In addition, we show that a rough fit to the integrated density of 
states with a single exponential: Nstates ~ exp(Ai5/0), where AE is the excitation energy, leads 
to a parameter, 0, exhibiting a universal scaling a la Thomas-Fermi: 9 = Z'^ N-^'^h{N / Z) , where 
h is approximately linear near N/Z — 1. 

PACS numbers: 32.30.-r, 32.10.Hq, 31.15.-p 



Energy levels of atoms and ions have been measured 
since the earliest times of Quantum Mechanics. At 
present, there is a huge amount of very precise measure- 
ments, which can be found, for example, in the critical 
compilation by the NIST^^. In these big volumes of data, 
sometimes universal relations remain hidden. In order to 
find them, you need a simple theory in which universality 
or scaling naturally appears. 

Thomas-Fermi theory^ was the first mean-field theory 
of atoms, the predecessor of more elaborated approaches, 
such as Density Functional TheorjSl It is based on a 
simple estimation of the kinetic energy, with the help of 
Pauli exclusion principle, and the introduction of self- 
consistency with the Poisson equation in order to take 
account of Coulomb interactions. Thomas- Fermi the- 
ory has proven to be a valuable tool for the qualita- 
tive and semi-quantitative understanding of rea l^* and 
artificial atoms. 

A very interesting aspect of Thomas-Fermi theory is 
the scaled form of physical magnitudes, which are very 
often respected by the exact quantum mechanical magni- 
tudes. This aspect has not been completely exploited so 
far. Recently, we showed, for example, that the density of 
low-lying excited states of artificial (quantum dot) atoms 
can be parametrized in a universal (scaled) way^, sug- 
gested by Thomas- Fermi theory. No similar parametriza- 
tion for neutral atoms or ions is available yet. 

In the present letter, we focus on the energy spectra 
of atomic ions, and show that, in some sense, the excita- 
tion spectrum is "universal", that is, may be described 
by scaled functions. In order to support the theoretical 
statements, we make extensive use of the NIST detailed 
compilation on atomic energy levels. 

Let us start considering the first ionization potential 
of ions, which sets the energy scale for the true discrete 
spectrum. In atomic units, where atomic ions are charac- 
terized by only two free parameters: the nuclear charge, 
Z, and the number of electrons, N, the ground-state en- 
ergy in the Thomas-Fermi approximation may be shown 
to satisfy the following scaled relation^: 
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FIG. 1: (Color online) Scaling in the ionization potential of 
ions with more than ten electrons. The curve {N/Z)~^^^ g, 
where g is given by Eq. (|3|, is also drawn (dashed line). 



Eg,{N,Z) = Z^N'/^'fiN/Z), 



(1) 



where / is a universal function of the variable N/Z . The 
ionization potential can be obtained from Eq. ([T|). In- 
deed, 



E,omz{N,Z) = Egs{N-l,Z)-Egs{N,Z) 

d_ 

~dN' 

« Z^N~^/^g{N/Z), (2) 



-rEgsiN,Z) 



where g is also a universal function of its argument. It 
may be easily verified for free electrons in a Coulomb 
field that, in the N ^ oo, N/Z limit, .g(0) = 
(3/2)1/3/3 «o.3816.fal 

The law given by Eq. ^ does not hold for neutral 
atoms. The ionization potential of neutral atoms, apart 
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Ar-like ions (N = 18) 
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FIG. 2: (Color online) The excitation spectra of Ar-like ions 
showing the Z'^ dependence of intra-shell excitations, and the 
exponential increase of Nutates with AE. 



FIG. 3: (Color online) Scaling of the temperature parameter, 
B. Ions with A*' > 10 and at least 100 measured lines are 
included in the plot. 



from shell effects, has a very smooth dependence on Z, 
and is close to that of Hydrogen, 0.5, in accordance to 
the fact that at long distances every singly ionized atom 
looks like a proton. 

However, the scaling law, Eq. is quickly reestab- 
lished in the cationic domain, Z > N. We show in Fig. 
[TJthe ionization potential of more than 380 atomic ions, 
taken from the NIST compilatioiP'. All of the systems 
with iV > 10 are considered, for which the scaled be- 
haviour is apparent. The nuclear charge, Z, spans the 
range from 12 (Mg) to 80 (Hg), although, as it can be 
seen from the table of available level^^, most of the data 
correspond to the interval 12 < Z < 37. 

The function g(x) should go to zero in the a; — s- 1 
limit, in accordance to the fact that, in the leading ap- 
proximation, the anionic instability border is located at 
N/Z ~ l.QJJ We fitted the experimental data with a 
two-point Pade approximantP^ interpolating between the 
a; — > and a; ^ 1 limits: 



(po +Pi.t)(1 ~ x) 
l + qix 



(3) 



where po — .g(0), and the fitting parameters are: pi = 
2.6555, qi — 11.4399. The fit gives a very good qual- 
itative description of the data. The three points at 
N/Z « 0.5, lying a little farther, correspond to the 
Xe XXVII, Xe XXVIII, and Xe XXIX ions, for which 
it seems that there is an overestimation in the reported 
valued. 

We stress that the fit given by the function g captures 
the general behaviour of Eioniz as a function of N and 
Z. Effects not included in Thomas-Fermi theory, such 
as shell or sub-shell filling, spin-orbit interactions, etc, 
could lead to relatively large (of the order of 10% or even 
larger) deviations from the universal behaviour, as seen 
in Fig. [1] 



Next , we attempt a global universal description of the 
excitation spectrum. Excitation energies, AS, cover the 
range {Q, Eioniz)- In the spectrum, we shall distinguish 
between intra-shell and inter-shell excitations, which ex- 
hibit different behaviours. 

In quality of example, we show in Fig. |2] the Ar-like 
series, from the Ar I atom {N = Z = 18) to the Mn VIII 
ion {N — 18, Z ~ 25). The y axis represents the total 
(accumulative) number of states for excitation energies 
below AE. The degeneracy of multiplets is explicitly 
taken into account. In the x axis, the excitation energy 
is scaled by Z^. 

First, we notice that from Ar I to Ca III the spectrum 
is restructured. This corresponds to reestablishing the 
Z — > oo sequence of Hartree-Fock orbitals: Is, 2s, 2p, 3s, 
3p, 3d,. . . , in place of the neutral atom sequenceHI. 

From Ca III on, it is evident that the low-lying spec- 
trum goes to a definite limit as Z increases. These are 
intra-shell excitations (in the = 18 system, the n = 3 
shell is open). This means that intra-shell excitation en- 
ergies scale as AE ~ Z'^, independently of the degree of 
ionization, N/Z. We notice that from Z oo perturba- 
tion theory one would expect intra-shell excitations to be 
of the order of Z, thus the observed behaviour is stress- 
ing that the lowest excitations have a non-perturbative 
character. 

Larger excitation energies correspond to inter-shell 
transitions, which are of the same order of Eioniz and, 
thus, exhibit the characteristic dependence on Z, N, and 
N/Z. 

Let us stress that for both the intra- and the inter- 
shell regimes, locally the number of states increases ex- 
ponentially with energy differences. This behaviour is 
quite generaP. It was first observed in nuclear sys- 
tems, where it is known as the "constant temperature 
approximation'''^^. In Fig. |2] dashed lines represent 
fits to the intra- and inter-shell excitations in the Mn 
VIII ion. The discontinuity of the slopes signals that the 
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mechanisms of formation of these states are differentia. 

A rough characterization of the spectrum can be given 
in terms of a single exponential: 

Nutates -exp{AE/e), (4) 

where we expect the temperature parameter, Q ^ Z^. 
The dependence of O on can be extracted from the 
Z ^ oo asymptote. In this limit and for closed shell 
ions, the gap and the number of electrons are expressed 
in terms of the principal quantum number of the last 
occupied shell as: AEi sa Z"^ jr? , and N w 2n^/3, re- 
spectively. The available number of levels at energies 
A£;i is Nutates - [2(n+ / {[2n'']\ [2(2n + 1)]!} ^ 
exp(an + bnlnn) ^ exp(cn), where a, b, and c are nu- 
merical coefficients. The last expression comes from the 
fact that Inn is a slowly varying function. Comparison 
with Eq. Q leads to 6 - Z^iV-^/s. 

In the spirit of Thomas-Fermi scaling, we assume for 
Q the relation: 

e Z^N-^^^h{N/Z), (5) 

where /i is a universal function. We fit the excitation 
spectra of ions with > 10 with a function given by Eq. 
Q. We study all of the ions with at least 100 measured 
lines, and use the first 100 lines to perform the fit. The 
results are shown in Fig. [3] 

Although Eq. Q gives a rough global characterization 
(the exponential dependence holds locally, with different 
slopes in different sectors of the spectrum), and there 
could be experimental problems such as missing lines, etc, 
the data shows scaling, with a function h approximately 
linear: 



0.045, with strong fluctuations related to shell effects. 
The reason for such a behaviour is also the proximity of 
the anionic instability border. 

In conclusion, we verified that Thomas-Fermi theory 
offers a fresh view to the energy spectra of atomic ions. 
The energy spectrum can be described in a "universal" 
way in terms of scaling functions. One of these func- 
tions describes the first ionization potential, which is the 
natural scale for the excitation spectrum. On the other 
hand, if the excitation spectrum is fitted by a single expo- 
nential, the temperature parameter is shown to approxi- 
mately exhibit also a scaling relation a la Thomas-Fermi. 

For neutral atoms, the scaling is broken because of 
the proximity of the anionic instability border. We shall 
further study this border seeking for positive uses of its 
presence. 

There are still many interesting points requiring fur- 
ther investigation. Intra- and inter-shell excitations 
could be more detailed described by means of constant- 
temperature approximations. Alternative characteriza- 
tions of the spectrum could also be employed. We have 
in mind, for example, the first moment, which is a mag- 
nitude with dimensions of energy. The energy region 
closer to the ionization potential could also be stud- 
ied. In this region, we expect the dependence Nstates ~ 
{Eiordz — AE)^^/'^ , coming from the existence of Rydberg 
states. A very interesting possibility is to write down an 
interpolation formula to spectroscopic accuracy for the 
ionization potential, which takes our g function as the 
starting point, but includes also shell-filling effects, spin- 
orbit interactions, etc. Research along these lines is in 
progress. 
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/i « 0.165 (1 - A^/Z). (6) 

Due to the reasons mentioned above, the observed dis- 
persion of points is natural. 

Similarly to the ionization potential, the parameter O 
for neutral atoms takes in mean a constant value, O « 
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